Effect of the on-site Coulomb repulsion on superconductivity 
in the boson-fermion model 
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We study the influence of the repulsive Coulomb interactions on thermodynamic properties of 
the boson fermion model with an anisotropic (d-wave, and extended s-wave) order parameter. Su- 
perconductivity is induced in this model from the anisotropic charge exchange interaction between 
the conduction band fermions (electrons or holes) and the immobile hard-core bosons (the localized 
electron pairs). The on-site Coulomb repulsion competes with this pairing interaction and hence is 
expected to have a detrimental influence on superconductivity. We analyze this effect in some detail 
considering the two opposite limits of: the weak and strong repulsion. A possible crossover between 
both these regimes is also discussed. 
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INTRODUCTION 

The boson fermion (BF) model describes a mixture of 
the narrow band fermions coupled to a system of the com- 
posite hard-core bosons. Initially, this type of an effective 
Hamiltonian has been invented for a system of itinerant 
electrons interacting with the local lattice deformations 
in the crossover regime, between adiabatic and antiadia- 
batic limits Later, the same model has been indepen- 
dently considered by a number of authors H, ||, ^, || as 
a possible scenario for a mechanism of the high tempera- 
ture superconductivity (HTSC). There are also attempts 
to apply a similar BF model to explain certain aspects of 
the Bose condensed atoms of the alkali metals ^ . 

This model reveals a rich physics both in its: nor- 
mal phase and the broken symmetry superconduct- 
ing/superfluid state. As shown in the mean field type 
studies 1^ there is a characteristic temperature bel- 
low which fermions are driven to the superconducting 
phase and simultaneously bosons start to Bose-condense. 
This result has been confirmed (neglecting the hardcore 
nature of bosons) by means of the shielded potential ap- 
proximation [0 and with a help of the renormalization 
group approach Moreover, when approaching the 
critical temperature from above, the pairing-wise corre- 
lations start to manifest themselves strongly. In partic- 
ular, they may give rise to a formation of the pseudogap 
in the fermion spectrum. This effect, known experimen- 
tally from a variety of measurements (see e.g. review pa- 
per 1^), provides a firm argument for application of this 
model to describe the HTSC materials. 

Pseudogap formation and its variation with a lowering 
temperature has been carefully investigated for the BF 
model using: a) the selfconsistent perturbative treatment 
of the boson-fermion coupling |ll| , b) the perturba- 
tive treatment of the kinetic processes (a la Hubbard I for 
the fermion hopping) with respect to the exact solution 
of this model in its atomic limit jl2|, c) the dynamical 
mean field theory (DMFT) equations which have been 
selfconsistently solved within the noncrossing approxima- 



tion for the auxiliary impurity problem p^ , d) and the 
renormalization group technique |^ . 

Many experimental data, especially the angle resolved 
photoemission spectroscopy Q, seem to suggest the 
anisotropic d-wave type structure of both: the pseudo- 
gap and true superconducting gap. However, there are 
also known some measurements, for instance the c-axis 
Josephson tunneling and the photoemission spec- 
troscopy on Bi2Sr2CaCu208-(-5 which provide the 
strong arguments for a nonzero s-wave ingredient of the 
order parameter. In a most realistic situation one can 
expect that the order parameter of the HTSC cuprates 
acquires a mixed s -I- d or s-\- id symmetry. Possibility 
for an appearance of the mixed symmetry superconduct- 
ing phase has been theoretically explored on quite the 
general grounds using a 2-dimensional electron system 
with the anisotropic potential V = Vs + Vd the arbi- 
trary (from weak to strong) attraction strength . So 
far, most of the studies of superconductivity within the 
BF model have been performed for the isotropic pairing 
interaction. Some attempts to analyze the d-pairing su- 
perconductivity together with a microscopic justification 
for introducing the BF type Hamiltonian can be found 
in the paper by Geshkenbein et al Very recently, a 
more formal way has been explored by Micnas et al . 

In this paper we shall investigate various kinds of the 
superconducting phase induced by an anisotropic poten- 
tial of the BF model in a presence of the Coulomb inter- 
actions between fermions. For simplicity we shall concen- 
trate only on a case of the on-site repulsion [/ > 0, where 
U = (m|-^|m) in the Wannier representation. In general, 
one expects that the on-site repulsion (which prevents 
fermions from forming the local pairs) would compete 
with the correlations induced by the boson-fermion cou- 
pling (this interaction is a driving force for the pairing in 
the BF model and is responsible for inducing the pseu- 
dogap at temperatures T* > T > Tc and the true su- 
perconducting gap when T < Tc). We shall address the 
following question: what is an extent of a detrimental 
influence of U on the pairing correlations ? 
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The abovementioned competition has been aheady 
studied in a normal phase of the BF model using the 
nonperturbative approach of the DMFT ||2^ . In our pa- 
per we shall investigate the anisotropic superconducting 
phase. To make our study feasible we assume that the 
pairing interactions are relatively weak (a meaning of this 
assumption is explained in the next section) and we try to 
estimate the influence of the Coulomb repulsion varying 
its intensity from the weak to strong interaction limits. 



THE MODEL 

Hamiltonian of the system under consideration consists 
of the two parts H = H^^ + H'^°"K First of them refers 
to the standard BF model Hamiltonian j5| 

= ^(ek-M)c^Ck,. + (As-2^)^6|6. 

k,cr i 

+ (^^k,q^>]jC„k+q/2,iCk+q/2,T + ^^-C-) (1) 

k,q 

and the second part denotes the on-site interaction be- 
tween fermions iJ^"™' = U^^niirii-^ . We use here the 
standard notations for annihilation (creation) operators 
of fermion Ci^^ (cJo-) with spin a and for the hard-core 
boson bi (bl) at site i of the 2 dimensional square lattice. 
The indices k and q in (|l|) denote the coordinates of the 
momentum space. We assume the tight binding disper- 
sion for fermions Ek = —2t{coskx + cos ky), and set the 
bandwidth D = 8t a.s a unit {D = 1). 

It is important to remark now that we let the boson- 
fermion exchange potential Wk.q to be anisotropic. As 
explained in the ref. ||l^ the low energy physics of this 
model is prevailed by bosons of small momenta |q| ~ 0. 
Usually, the magnitudes of the superconducting gap in 
HTSC materials are of the order of several meV (which 
is ~ lO^'^ of the bandwidth D). It is thus reasonable 
to assume that the pairing potential Wk,q, which estab- 
lishes the energy scale for Tc and Asc{T = 0), is small as 
compared to D. In such a case the following mean field 
decoupling is justified 

^ Wk,q&qC-k+q/2aCk+q/2,T - 

q 

Wk,q=0 ((6q=0>* C-kjCk^T + &q=0 (C-k,iCk,T)) • (2) 

We further write down the anisotropic potential Wk.o as 
a product 1 19 

^'k.o = g <^k, (3) 

where g characterizes the interaction strength and 0k 
stands for the dimensionless factor which has to reflect 
the fourfold symmetry of Cu02 planes of the HTSC 



cuprates. In general the anisotropy factor (f>\^ can be rep- 
resented as 

<f>k = ctf) + as (cos + cos ky) + ad (cos k^ — cos ky) (4) 

and the coefficients ao.s.d denote a relative contribution 
of the isotropic, the extended s-wave, and the d-wave 
part into the order parameter of the superconducting 
phase. They should be adjusted depending on a spe- 
cific material. If for example ao ^ and ^ we 
would have the order parameter of a mixed s + d, or 
s+id symmetry. Since our main interest is focused on 
the competition between the Coulomb interaction and 
the superconductivity we further consider for a clarity 
only the pure extended s- or the rf-wave symmetries when 
0k — cos kx ± cos ky . 

After the mean field decoupling for the boson-fermion 
interaction (|^) we are left with an effective Hamiltonian 
composed of the separated fermion and boson contribu- 
tions I! H ~ iJ^ + iJ^ 

= ^-^kCk^^Ck,^ + C/^nj,xni,t 

k,(T i 

+ P </'kc[^cLk,x + h-c.) , (5) 



H"" = J2(Eoblh + gxbl+gx*b, 



(6) 



We introduced here the abbreviations for energies mea- 
sured from the chemical potential $k = £k ^ 1^-, Eg = 
As — 2/i and for the two order parameters x — 
Ek'/'k (c-k,iCk,T), P= (&q=o)/\/A^= {bt). 

We can easily solve the hard-core boson part of the 
problem. For a given site i one finds the true eigenstates 
using the unitary transformation 

\A>i = cos(a)|0 >j +sin(a)|l >j (7) 
\B>, = -sin(a)|0 >, +cos(a)|l >i (8) 

such that tan (2q;) — {~2gx)/ Eq, where |0 >i, |1 >i refer 
correspondingly to the empty and to the singly occupied 
(by the hard-core boson) site z. In a straightforward cal- 
culation we can determine the expectation values for the 

number operator = (blbi^ and for the order pa- 
rameter p j|, [l9| 



l_Eo 
2 47 
gx 

— '■ — tanh 
27 



tanh 



knT 



ksT 



(9) 
(10) 



where 7 = \/ {Eo/2)^ + \gx\^ and ks is the Boltzmann 
constant. 
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WEAK INTERACTION LIMIT 

First, we consider the weak couphng hmit when 
U is fairly smaller than D. We are in a position 
to utilize then the Hartree Fock Gorkov lineariza- 
tion for the on-site interaction ni^|n.i_| ~ n^^j {ni^^) + 

(n^,x)n,,| + c^^id^^^ c\^c\ ,^ (cj,iC»,T)- Hamil- 

tonian of the fermion subsystem (^) reduces then sim- 
ply to the BCS structure ~ ^ ^kcj^. „c\^,a + 

Ek (^k -^^^^jcLka + ^-c-) ^it^ = a + t^^^'/a (we 

assume a paramagnetic state = ('T'i,!) = A 

role of the effective gap parameter is played here by 



Ao + g p <^k, 



(11) 



where the isotropic part is given by Aq = J7 (ci^|Ci_|). 
Standard methods of the solid state theory give the fol- 
lowing equations for expectation values 



(c-k,iCk,T; 



1 — -|r- tanh 
k 

■ tanh 



-Al 
2Ev 



£^k 



(12) 



(13) 



Ajf^-^'P in 



with a typical gaped spectrum = 
the superconducting phase. 

It is worth mentioning that in a case of the d-pairing 
(i.e. for 0k = cosfcj; — cosky = ?7k) the isotropic compo- 
nent Aq of the gap parameter (|l^) does identically van- 
ish. To prove this let us substitute (|l^) into the definition 
of Aq = Uj2k ('^-k4Ck,t) to obtain 

^ 2Ei, \2kBT) ^ ' 



Aq: 



Since integration over the Brillouin zone of the part con- 
taining 77k gives zero so, for [/ > 0, equation ([ij) has 
the only possible solution Ag = 0. It is not surprising 
because the repulsive interactions by themselves are not 
able to induce the on-site fermion pairs. 

If the boson fermion potential (^ is isotropic or takes 
a form of the s-wave ((/ik = cos kx + cos ky) then in gen- 
eral Aq 7^ 0. From a consideration similar to the one 
discussed above (equation ( p^ ) is valid except that 77k 
should be replaced by 0k) we can determine a relative 
ratio Ao/gp. The extended s-wave gap parameter is now 
given by 



A. 



9 P 



{eff) 



Ek jUM'^Ek) tanh {E^/2kBT) 
1 + Ek {UI2E^) tanh {E^/2kBT) 



(.15) 



For the isotropic boson fermion potential (0k = 1) the 
equation (nS) simplifies further to give a k-independent 



gap A(^//) = 5P [1 + Ek iU/2E^) tanh {E^/2kBT)]-^ < 
gp. This expression explicitly shows a detrimental role 
of the on-site repulsion on the isotropic superconducting 
phase. Such a problem has been previously addressed 
[0, ^ neglecting the hard core nature of bosons and 
using the RPA treatment for the Coulomb repulsion. 



STRONG INTERACTION LIMIT 

In a case of the strong interactions {U > D) we make 
use of the slave boson technique proposed by Kotliar and 
Ruckenstein |22). For simplicity we shall consider here 
only the extreme limit [/ — > 00. 



We represent the fermion operators as 



4k 



and = fja^i' where the auxiliary boson operator 
Ui (a|) refers to the annihilation (creation) of the empty 
state at site i and fermion operator fi^cr ifja) corresponds 
to annihilation (creation) of the singly occupied site i 
with spin a. No double occupancy is allowed and this 
can be formally expressed via the local constraint a^a^ -I- 

Ecr /j.cr/i,"" ~ 1- 

Using the real space (Wannier states) representation 
we can rewrite the Hamiltonian (|^) in terms of the new 
operators as 



(16) 



We used here the identity aja| = 1 and the last 
term takes account of the local constraint (Ai stands for 
the Lagrange multiplier). Vi,j is the exchange potential 
whose Fourier transform is given by (^) and, as usually, 
ti,j denotes the hopping integral. 

Next, we approximate ( p^ by: (i) replacing the slave 
boson operators by their expectation values which are 
assumed to be site independent Oi ~ (a^) ~ r, and (ii) 
replacing the local multipliers by the global one Xi ~ A. 
In this (mean field) approximation for the slave bosons 
one obtains 



k,cr 



(17) 



The global parameters A, r are determined from a min- 
imization of the total energy (H). This criterion leads 
to 



1 



(18) 
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E 



2Relgp*y6^ (/-k,i/k,T) 



k 



(19) 



As can be noticed from the equation (|T^) Hamiltonian 
of the fermion subsystem is again reduced to the 
BCS structure. We thus have the same solution for the 
expectation values as given in ( p^ , p^ with a difference 
that now 



6 
Ak 



r ei^- n + X , 



Both, the effective fermion bandwidth D^'^^^^ 



(20) 
(21) 

r^D 

and the effective pairing potential V^'^-'^-^-' = r^94>'kj re- 
duce down to zero when fermion occupation approaches 
the half-filhng. Under such circumstances the system is 
driven into the Mott insulating state. 



THE CROSSOVER 

Finally we consider here a regime of the intermediate U 
for which we adopt the procedure used earlier by us |^ in 
a context of the extended Hubbard model. We introduce 
the Nambu representation = (cj^ -|.c_k,j,), ^k = (^k)^ 
and, as a first step, determine the unperturbed Green's 

fimction G°(k,uj) = ^^^'k;*k^^ neglecting the inter- 
action U in the fermion Hamiltonian (|^) 



i^^n - £.k - gp*<j)k 

-gp4>k ii^n - S,k 



(22) 



Next, we compute the dressed Green's function using the 
matrix Dyson equation 



G """(k, ZWn 



Gq ^ (k, iuj„ 



(23) 



In order to proceed we simplify the selfenergy matrix by 
the following Ansatz [Esl 



E(k, iuj„ 



EAr(k,itj„) U (^cl^cl^'^ 



(24) 



Without specifying the diagonal elements of ( p^ we de- 
note Eii(k, iujn) by I]Ar(k, iLUn) and make use of the iden- 
tity E22(k, iwn) = — Sii(k, — iw„). The off-diagonal ele- 
ments are approximated by us by a contribution corre- 
sponding to the result deduced from the mean field type 
theory discussed in the section HI. Channel of the super- 
conducting correlations is treated by us in ( p^ approxi- 
mately. 

If one knew I]Ar(k, icjn) then the needed expectation 
values can be found according to the standard field the- 
oretical relation (AB) = X^^oo ((^' ' where 



/? = l/ksT. In particular, we obtain [p3[ 



4^CkT 



(c-kiCkt; 



i^^n + ^k + ^N{-i^n) 



^ Kc^„-ek-SA'(ic^„)P-MA(,^^/)|2 

(25) 



Aeff) 



ZC^„-ek-SA'(*C^„)P + |AL'-^^^|2' 



(26) 



where again A^'^'^''^'' is given by (^jj). 

It is worth noticing that I]jv(k, iuJn) has a meaning of 
the normal phase selfenergy for the standard Hubbard 
model. Of course there is no exact solution for S n avail- 
able so far except maybe from the numerical exact diag- 
onalization or the Quantum Monte Carlo studies. How- 
ever, depending on a magnitude of the on-site interac- 
tion, one can use various approximate estimations. Let 
us point out the few possibilities. 

Starting from the weak interaction limit the simplest 
substitution for E^v is the mean field value Un^ /2 as 
discussed in section HI. With an increase of U one can 
proceed by including some higher order corrections, like 
for example of the second order in U p^ . Going toward 
the Mott transition regime U = Ucr ~ D (for the half- 
filled fermion system) one could work for example with 
the so called alloy analogy approximation (AAA) ||2^ . In 
a more subtle way one could estimate the momentum in- 
dependent selfenergy Yi]^{lo) by adopting the dynamical 
mean field theory (which becomes exact in the limit of 
infinite spatial dimensions). The strong interaction case 
can be described in a satisfactory way either with AAA, 
DMFT or the simple Hubbard I approximation. Here 
we apply the AAA procedure to study a regime near the 
Mott transition. 



NUMERICAL RESULTS 

As far as the Coulomb repulsion is concerned we expect 
that its effectiveness should strongly depend on concen- 
tration of fermions. In a regime of small fermion con- 
centration (dilute limit) this interaction should not be 
very efficient and, in particular, it would not affect the 
superconducting type correlations induced by the boson 
fermion exchange. On the other hand, we expect that the 
strongest effects of the on-site Coulomb repulsion might 
appear for a system with a nearly half filled fermion sys- 
tem nF = 1. 

In absence of the Coulomb interactions, the super- 
conducting phase (isotropic or anisotropic one) of the 
BF model is formed when the fermion concentration is 
properly adjusted. The Fermi energy ep must be close 
enough to the boson level because only then the charge 
exchange between the hard-core bosons and fermion pairs 
can induce the long range coherence ||2^. So, the needed 
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Ab/2 



A3=-0.6 




FIG. 1: A schematic illustration of the two distinct critical 
fermion concentrations: (top) np — 0.4 when boson energy 
As/2 = -0.3, and (bottom) = 1 when Ab/2 = 0. The 
shaded areas show the occupied fermion states at zero tem- 
perature. 



concentration of fermions is roughly given by = 
2 /f^^^^^^ /5o(e)rfe, where po{£) denotes the density of 
states of the free (nonintcracting) fermion system. 

In this section we present results obtained numerically 
for the BF model in the two distinct cases when the crit- 
ical fermion concentration is: a) small, which takes place 
when Ab/2 lies fair aside the center of fermion band, b) 
close to half- filling, when Ab/2 is located in the center of 
fermion band. We thus choose the two following values 
Ab/2 = -0.3 and 0, see figure (0). 

We take the boson fermion potential 5 = 0.1 in all 
the results discussed below. Figure (||) shows the critical 
temperature T,. of the d- wave superconducting phase for 
several values of U . This type of superconductivity is en- 
hanced near the half-filled fermion system similarly as in 
the extended Hubbard model (25j. In agreement to our 
expectations, the Coulomb repulsion only weakly reduces 
Tc in a dilute fermion system. However, for Ab = we 
notice a considerable reduction of Tc or even a disappear- 
ance of superconductivity for np — 1 (total concentration 
is then ntot — np + 2nB = 2) when U exceeds the critical 
Mott transition value Ucr- For such a strong potential U 
the d-wave superconducting phase is restricted to a nar- 
rower regime of the total concentration, such that there 
are no doubly occupied fermion states on a given lattice 
site (remember we are considering the intersite Cooper 
pairs). 

Figure (^ illustrates the effect of U on the extended 
s-wave phase. In a dilute regime we notice almost iden- 
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FIG. 2: Transition temperature Tc into the d-wave super- 
conducting phase as a function of the total concentration ntot 
for the two representative boson level values As = —0.6 (top) 
and As = (bottom). Curves corresponding XoU = 0.1 and 
0.2 where obtained from the mean field approximation, while 
17 = CO from the salve boson study. 



tical values of Tc for both the d- and extended s-wave 
superconducting phases. Also the influence of U is there 
very similar. A remarkable difference appears for Ab = 
when the fermion system is near the half-filling np — \. 
Critical temperature T^ of the s-wave phase is then 3 
times smaller as compared to the c?-wave phase. The sys- 
tem is thus less susceptible for the s type pairing near 
np = \ (for Ab = it corresponds to ntot — 2). 

One notices also some "peculiar" behavior oiTcijitot) 
for the extended s-wave phase cased by an increasing 
strength of I] . With a small increase of U the whole dia- 
gram is somewhat shifted and simultaneously the optimal 
value of Tc slightly increases. This overall shift is caused 
by the Hartree term Unp/2 (see the section III) which 
effectively pulls up the fermion band with respect to the 
boson energy level. By comparing the curves correspond- 
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0.01 




0.01 




FIG. 4: Transition temperature Tc of the d- and the extended 
s-wave phases for the half-filled fermion system riF = 1- Both 
superconducting phases disappear when the Mott transition 
is approached. The inset shows the density of states at the 
Fermi energy for a normal phase obtained by the Alloy Anal- 
ogy Approximation. 



on the dynamical mean field theory usually yields Ucr 
larger than 1 Nevertheless, a qualitative behavior 

presented in figure remains valid. As we see by com- 
paring to the figures the AAA treatment properly 
interpolates between U ~ Q and U = oo limits. 



FIG. 3: Transition temperature Tc into the extended s-wave 
superconducting phase for several values of the on-site repul- 
sion, as indicated. The top figure corresponds to As = —0.6 
and the bottom one to A_g = 0. 



ing to [/ = in the upper and bottom panels of figure (j^) 
we realize that such shift is responsible for enhancing the 
s-wave type superconductivity, but only when U is safely 
smaller than D. A further increase of the Coulomb inter- 
action U proves to be detrimental on superconductivity 
(independently of a symmetry of the order parameter) 
because the fermion subsystem is driven into the Mott 
insulating state. 

To analyze in more detail a pronounced effect of the 
Coulomb interaction on superconductivity for the half- 
filled fermion system riF = 1 {ntot = 2) we show in figure 
the dependence of Tc on U. The results have been 
obtained by means of the Alloy Analogy Approximation 
mentioned in section III and discussed earlier by the same 
author in Ref . ||2^ . For the dim — 2 tight binding disper- 
sion we determine the Mott transition at Ucr — 0.54 in 
units of the initial fermion band. This value is probably 
underestimated. The most credible determination based 



CONCLUSIONS 

Summarizing, we have investigated the anisotropic su- 
perconductivity within the boson fermion model in a 
presence of the Coulomb repulsion between fermions. 

(a) In a dilute regime of the fermion concentration the 
effect of the Coulomb repulsion proves to be rather weak. 
For both the d- and extended s-wave phase we observe up 
to 25 % reduction of the optimal Tc value when U — s- oo. 
Both anisotropic phases survive, even in the limit of in- 
finitely strong Coulomb repulsions. 

(b) In the nearly half-filled fermion system we observe an 
enhancement of the c?-wave superconducting phase and a 
simultaneous suppression of the s-wave phase until the 
interaction U is small. 

(c) Around the Mott transition Ucr — 0.54 both phases 
are reduced to the concentration regime np < 2nB < 
1. Still, the superconductivity is able to survive at suffi- 
ciently large hole concentration h = 1 — np > 0. Such a 
case is relevant for a description of the HTSC materials 
and boson fermion model seems to capable to reproduce 
qualitatively the phase diagrams known for these mate- 
rials. 

Among the problems which are not addressed in this 
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paper there is a very intriguing question: what happens 
to pseudogap of the normal phase (discussed earher in 
the Refs |, ^ |l^, |l| and in ^ in a presence of 
Coulomb interactions ? Consideration of this subject is 
in a progress and the results shall be discussed elsewhere. 

Author kindly acknowledges helpful discussions with 
J. Ranninger and K.I. Wysokihski. Partial support has 
been provided by the Polish Committee of Scientific Re- 
search under project No. 2P03B 106 18. 
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